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Introduction 

The  ability  of  near-infrared  (NIR)  light-based  techniques  to  noninvasively  image 
and  analyze  tissue  structure  and  function  promises  their  great  potential  for  detection 
and  diagnosis  of  breast  cancer.  Optical  diagnostic  techniques  allow  us  to  not  only 
enhance  the  existing  capabilities,  but  to  eliminate  the  need  for  physical  biopsies.  In 
addition,  optical  imaging  is  inexpensive  and  portable,  which  indicates  that  optical 
imaging  could  be  an  ideal  candidate  for  routine  breast  screening.  However,  since  the 
scattering  properties  of  tissues  convolute  re-emitted  NIR  signals,  the  extraction  of 
pertinent  information  continues  to  remain  elusive.  An  understanding  of  light  propagation 
and  light-tissue  interaction  is  required  before  the  optical  technologies  can  substantially 
impact  diagnostic  medicine. 

Research  efforts  in  the  Biomedical  Optics  Laboratory  at  Clemson  University  are 
focused  on  the  biophysics  of  light  propagation  and  light-tissue  interaction  in  order  to 
engineer  appropriate  approaches  for  noninvasive  breast  imaging  and  spectroscopy. 
Specifically,  we  are  developing  indirect  optical/fluorescence  approaches  using  photon 
migration  measurements  in  the  continuous-wave  and  frequency  domains.  These 
indirect  optical/fluorescence  approaches  or  image  reconstructions  are  computationally 
based  on  the  powerful  finite  element  methods.  A  CCD-based  optical  spectroscopic 
imaging  system  is  already  operational  in  our  laboratory  for  continuous-wave 
tomographic  photon  migration  measurements,  while  a  frequency-domain  system  has 
also  been  constructed.  Using  these  optical  systems  coupled  with  our  finite  element 
based  reconstruction  algorithms,  we  will  be  able  to  extract  spatial/spectroscopic  maps 
of  tissue  optical  properties,  lifetime  and/or  yield  of  endogenous  and  exogenous 
fluorescent  probes.  Since  metabolic  tissue  states  can  be  identified  by  our  approaches, 
diagnostic  information  is  also  obtained  in  addition  to  detection  of  tumor. 

This  Career  Development  application  for  support  of  Dr.  Huabei  Jiang  will  facilitate 
the  establishment/continuation  of  these  research  activities.  Interdisciplinary  interactions 
with  the  Greenville  Hospital  System  (Greenville,  SC)  will  be  enhanced,  which  insures 
the  direction  of  research  towards  a  clinically  pertinent  and  feasible  system. 

Body 


This  report  describes  work  accomplished  during  the  fourth  year  of  a  proposed 
four-year  study  (with  one-year  no-cost  extension).  This  Career  Award  supports  Dr. 
Jiang’s  research  on  optical  and  fluorescence  imaging  using  both  continuous-wave  and 
frequency-domain  measurements.  The  focus  of  the  proposed  work  in  Year  4  is  the 
implementation  and  evaluation  of  image  enhancing  schemes  in  the  areas  of  dual 
meshing  and  total  variation;  continued  phantom  studies  of  image  reconstructions  with 
dye-free  or  dye-laden  phantom  background;  continued  implementation  of  3D 
reconstruction  algorithms. 

Software  work:  We  have  implemented  and  evaluated  important  image  enhancing 
schemes  in  the  areas  of  dual  meshing  and  total  variation.  Extensive  simulation  and 
phantom  experiments  have  been  conducted  to  complete  these  software  developments. 
We  have  also  started  to  implement  the  3D  reconstruction  algorithms  based  on  these 
image  enhancing  schemes. 
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Phantom  Experiments:  Using  our  frequency-domain  imaging  system,  we  have 
conducted  extensive  phantom  experiments  for  both  optical  and  fluorescence 
reconstruction  with  dye-free  background  and  dye-laden  background.  Our  experimental 
setup  used  was  the  automated  multi-channel  frequency-domain  system  mentioned 
above.  The  system  employed  a  radio-frequency  intensity-modulated  near-infrared 
beam.  The  laser  beam  was  sent  to  the  phantom  by  16  fiber  optic  bundles  coupled  with  a 
high  precision  moving  stage.  The  diffused  radiation  was  received  by  another  16  channel 
fiber  optic  bundles  and  delivered  to  a  thermo-electric  cooled  PMT.  A  second  PMT  was 
used  to  record  the  reference  signal.  These  PMTs  were  supplied  at  a  radio-frequency 
modulated  current  with  0.1-1  KHz  shift.  The  intermediary  frequency  signal  obtained  from 
the  PMTs  was  processed  using  a  National  Instruments  board.  For  every  source 
position,  16  measurements  for  each  detector  were  made,  taking  alternatively  100  ms 
samples  for  sample  and  reference  signals.  Fluorescence  signals  were  obtained  through 
an  830nm  or  690nm  interference  filter  placed  in  front  of  the  detection  PMT.  dc,  ac 
intensity  and  phase  shift  between  reference  and  sample  signals  were  obtained  using 
FFT  Labview  routines.  The  total  data  collection  time  for  256  measurements  was  8 
minutes. 

Solid  phantom  was  used  to  mimic  the  human  tissue.  It  was  made  of  agar, 
Intralipid,  black  ink  and  fluorescent  dyes.  The  absorption  and  the  reduced  scattering 
coefficients  are  linear  with  the  ink  and  Intralipid  concentrations,  respectively.  While  ICG 
and  DTTCI  dyes  were  continually  used  in  these  experiments,  we  have  also  used 
oxygen-sensitive  dyes  including  Sn(IV)Chlorin-e6-CI2-3Na  (SCCN).  Agar  was  used  to 
make  the  phantom  solid.  The  absorption  and  fluorescent  emission  peaks  of  this  oxygen- 
sensitive  dye  is  about  660  nm  and  690nm,  respectively.  The  solid  phantom  consisted  of 
a  cylindrical  background  and  a  cylindrical  heterogeneity.  These  experiments  show  that 
we  may  be  able  to  obtain  the  tissue  oxygen  concentration  images  by  the  use  of  oxygen- 
sensitive  dyes.  This  could  provide  optical  imaging  a  unique  capability  to  distinguish 
between  benign  and  malignant  breast  tumors  because  the  oxygen  concentration  in 
these  two  types  of  tumors  often  is  distinct.  Some  of  the  image  results  have  been 
published  in  a  peer-reviewed  journal  (see  the  manuscript  provided  in  the  Appendix  to 
this  Summary  Report). 

Key  Research  Accomplishments 

1.  We  have  implemented  and  evaluated  the  2D  image  enhancing  schemes  especially 
in  the  areas  of  dual  meshing  and  total  variation.  We  have  also  started  to  implement 
the  3D  reconstruction  algorithms  based  on  these  image  enhancing  schemes. 

2.  We  have  conducted  considerable  phantom  experiments  for  optical  and  fluorescence 
lifetime  imaging.  The  successful  experiments  have  confirmed  the  imaging  capability 
of  our  reconstruction  software. 

Reportable  Outcomes  (see  the  Appendix  to  this  Summary  Report) 
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Conclusions 

We  have  made  a  significant  progress  that  has  fulfilled  most  of  the  statement  of 
work  proposed  for  Year  4  of  this  project.  We  believe  that  we  will  be  able  to  fulfil  or 
exceed  the  work  statement  for  Year  4  within  the  extended  year. 


Appendix 

E.  Shives,  Y.  Xu,  and  H.  Jiang,  "Fluorescence  lifetime  tomography  of  turbid  media  based 
on  an  oxygen-sensitive  dye,"  Opt.  Express  10,  1557-1562  (2002). 
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Abstract:  We  present  for  the  first  time  experimental  images  of  fluorescence 
lifetime  distribution  using  model-based  reconstruction.  The  lifetime 
distribution  in  our  phantom  experiments  was  realized  through  using  an 
oxygen-sensitive  dye  [Sn(IV)Chlorin-eg-Cl2-3Na  (SCCN)]  whose  lifetime 
varied  with  the  oxygen  concentration  provided  in  the  target  and  background 
media.  The  fluorescence  tomographic  data  was  obtained  using  our  multi¬ 
channel  frequency-domain  system.  Spatial  maps  of  fluorescence  lifetime 
were  achieved  with  a  finite  element  based  reconstruction  algorithm. 

©2002  Optical  Society  of  America 

OCIS  codes:  (110.6960)  Tomography;  (170.3010)  Image  reconstruction  techniques;  (170.3650) 
Lifetime-based  sensing;  (170.6280)  Spectroscopy,  fluorescence  and  luminescence 
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1.  Introduction 

Recently  the  idea  of  reconstruction-based  fluorescence  diffusion  tomography  has  been 
developed  for  clinical  applications  such  as  breast  cancer  detection  and  tissue  functional 
mapping  [1-4].  This  new  imaging  approach  relies  on  the  fact  that  the  lifetime  of  fluorophores 
in  tissue  can  potentially  provide  tissue  functional  information  such  as  tissue  oxygenation,  pH, 
and  enzyme  [2-5],  It  is  also  based  on  the  fact  that  the  fluorophores  may  preferentially 
accumulate  in  tumors,  hence  providing  enhanced  sensitivity  for  cancer  detection  [6,7],  In  this 
type  of  indirect  imaging  method,  a  model-based  reconstruction  algorithm  is  crucial  which 
allows  for  the  formation  of  a  spatial  map  of  fluorophore  concentration  and/or  lifetime  in 
tissue. 

This  paper  explores  the  possibility  of  fluorescence  lifetime  tomography  (FLT)  based  on 
oxygen-sensitive  dyes.  Hohenberger  et  al.  [8]  have  recently  shown  that  the  oxygen 
concentration  for  healthy  and  cancerous  breast  tissue  is  significantly  different.  They  reported 
the  values  of  oxygen  concentration  for  healthy,  fibroadenomas  (benign  lesions)  and  malignant 
lesions  as  56,  63  and  44  mmHg,  respectively.  Fibroadenomas  have  higher  oxygen  content 
than  the  surrounding  tissue  because  the  solubility  for  oxygen  of  the  fluid  inside  is  presumably 
higher.  Conversely,  a  malignant  tumor  has  lower  oxygen  content  than  the  surrounding  tissue 
because  of  the  high  rate  of  mitosis.  Mitosis  requires  a  lot  of  energy  from  the  cell  that  raises  the 
cells  metabolic  rate.  As  the  metabolic  rate  of  the  cell  increases,  oxygen  use  goes  up  and 
therefore  the  oxygen  content  in  the  malignant  tumor  is  less  than  the  surrounding  tissue.  It  has 
also  been  shown  in  the  literature  [9-12]  that  the  lifetimes  of  some  fluorescent  dyes  are 
sensitive  to  the  amount  of  oxygen  present  in  solution/tissue.  The  more  oxygen  there  is,  the 
lower  the  lifetime:  oxygen  concentration  and  lifetime  is  quantitatively  correlated  by  the  well- 
known  Stem-Volmer  equation  [13].  Thus  tissue  oxygen  concentration  maps  can  be  derived 
from  lifetime  images  obtained  using  FLT. 

To  date  fluorescence  reconstructions  are  largely  limited  to  imaging  of  fluorophore 
concentration  [14],  or  fluorescence-enhanced  absorption  [15],  In  this  paper  we  present 
successful  quantitative  reconstruction  of  lifetime  distributions  from  frequency-domain 
measurements  using  tissue-like  phantoms  containing  SCCN  dye  with  varied  oxygen  content  in 
the  target  and  background.  To  the  best  of  our  knowledge,  the  results  shown  here  are  the  first 
experimentally  reconstructed  images  for  which  recovery  of  lifetime  profile  has  been  achieved 
with  absolute  ac  excitation  and  fluorescent  emission  data. 

2.  Methods  and  Materials 
2.1.  Reconstruction  Algorithm 

Based  on  a  set  of  coupled  diffusion  equations  that  describe  the  propagation  of  both  excitation 
and  fluorescent  emission  light  in  tissues,  our  reconstruction  approach  casts  image  formation  as 
a  nonlinear  optimization  problem  in  which  the  optimization  parameters  are  optical  and 
fluorescent  properties.  While  the  details  of  our  reconstruction  algorithm  can  be  found  in  Ref. 
4,  we  outline  a  special  formulation  for  the  phantom  study  reported  here.  In  frequency-domain, 
it  is  known  that  propagation  of  both  excitation  and  fluorescent  emission  light  in  tissues  can  be 
described  by  the  following  coupled  diffusion  equations  [3,16]: 
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(1) 


V-[Dx(r)VOx(r,co)]-  jxa  (r) - <J>x(r,a>)  = -S(r,co) 

c 


V  [D„(r)V®„(r,oi)]-  4>„(r,ffl) - <2> 

where  Oxm  is  the  photon  density  for  excitation  (subscript  x)  or  fluorescent  light  (subscript 
m),  Dx  m  is  the  diffusion  coefficient,  ^  is  the  absorption  coefficient  due  to  contributions 

from  both  non-fluorescing  chromophores  and  fluorescent  dye,  Hax_m  is  the  absorption 
coefficient  for  the  excitation  light  due  to  contribution  from  fluorescent  dye,  (0  is  the 
modulation  frequency,  c  is  the  velocity  of  light  in  the  medium,  andT|  and  T  are  the 
fluorescent  quantum  yield  and  lifetime,  respectively.  S(r,(0)  is  the  excitation  source  term  in 
(1)  which  for  a  point  source  can  be  written  as  S  =  S08(r-r0),  where  S0  is  the  source 
strength  and  8(r-  r0)  is  the  Dirac-delta  function  for  a  source  at  r0.  The  diffusion  coefficient 
can  be  written  as  Dx,m(r)  =  1/3[lJ'ax.m(r)  +  M'sjm(r)]  where  ^L,(r)  is  the  reduced 
scattering  coefficient.  In  this  study  we  use  the  non-zero  photon  density  boundary  conditions 
(BCs)  or  Type  III  BCs:  — E>xmVOx  m  •  n  =  0tOx  m,  where  n  is  the  unit  normal  vector  for 
the  boundary  surface  and  CL  is  a  coefficient  that  is  related  to  the  internal  reflection  at  the 
boundary. 

Making  use  of  finite  element  discretizations,  we  can  obtain  two  matrix  equations  for  Eqs. 
(1)  and  (2)  and  realize  other  derived  matrix  relations  through  differentiation,  which  lead  to  a 
set  of  equations  capable  of  inverse  problem  solution: 

[A,,  ]{«>,.„  }={b„}  (3) 

/3z}=  {»„  /ax}-PA,.„  /3xl{®,,}  (4) 

+  Xl)ix  =  3I„(®;,.  -«>:,)  (5) 


elements 


(ax,m)ij  =<  — Dx  mV\|/j  -  Vxj/j  —  ()_l ^ - )\|/j\j/i  >  where  <  >  indicates  integration 

0 

over  the  problem  domain;  \|/j  and  \)/j  are  locally  spatially-varying  Lagrangian  basis 

M 

functions;  the  entries  in  column  vectors  { bx>m }  are  (bx)j  =  — (S\|/;)  +  a£(Ox)i\|/j\l/ids 

j=i 

and  (bm)i  =-<nna  m  J(Ox)jVj¥.-^^->+aj;(a>ra)j|\lfj\|/ids  where 


and  (bm)i=-<TFa_Z(<l>x)j¥jVi7 

j=i  1 

<j”  expresses  integration  over  the  boundary  surf 


<j  expresses  integration  over  the  boundary  surface  where  Type  III  BCs  have  been  applied.  % 
expresses  Dx  ,  fla  ,  T(jla  ,  or  X ;  3  x  m  is  the  Jacobian  matrix  that  should  be  formed 
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from  30  x  m  /  3%  at  the  boundary  measurement  sites;  A%  is  the  update  vectors  for  the 
optical  and  fluorescent  property  profiles;  I  is  the  identity  matrix;  %  may  be  a  scalar  or  a 
diagonal  matrix;  and  ^x’m  =  •"(^Dn)1  and  (0°m)j  and 

(Oxm)j,  respectively,  are  observed  and  calculated  photon  density  for  i=l,  2,  M 


boundary  locations.  Note  that  to  estimate  Dx ,  (Ia>  ,  T|(i,a ^  ,  and  t  spatially  we  expand 

these  quantities  in  a  similar  fashion  to  O  as  a  finite  sum  of  unknown  coefficients  multiplied 
by  the  locally  defined  Lagrangian  basis  functions.  In  FLT,  images  are  formed  by  iteratively 
solving  Eqs.  (3)-(5)  and  updating  the  optical  and  fluorescent  property  distributions  from 
presumably  uniform  initial  estimates  of  these  properties. 


Fig.  1.  (a)  and  (b)  Experimental  geometries  under  study,  with  an  off-centered  target  7mm  away 
from  the  boundary  (a)  and  an  off-centered  target  9  mm  away  from  the  boundary  (b).  (c)  Gas 
delivery  system. 


2.2.  Experiments 

Our  experimental  setup  is  an  automated  multi-channel  frequency-domain  imaging  system  that 
has  been  described  in  detail  elsewhere  [17,18].  In  this  system,  an  intensity-modulated  635-nm 
light  from  a  30  mW  diode  laser  (Thorlabs)  is  sequentially  sent  to  the  phantom  by  sixteen  1- 
mm  fiber  optic  bundles.  For  each  source  position,  the  diffused  light  is  received  at  sixteen 
detector  positions  along  the  surface  of  the  cylindrical  phantom  and  sequentially  delivered  to  a 
photomultiplier  tube  (PMT)  (Hamamatsu  R928).  A  second  PMT  is  used  to  record  the 
reference  signal  (<5%  of  the  modulated  light  from  the  diode  laser).  The  two  PMTs  are 
supplied  by  the  radio-frequency  (RF)  generators  (IFR  Americas  2032A)  at  a  modulated 
current  of  100  MHz  with  1.0  KHz  shift.  Both  source  and  detector  fiber  bundles  are  firmly  held 
by  a  black-painted  metal  ring  structure.  The  multiplexing  of  the  source/detector  fibers  are 
accomplished  by  two  moving  stages  (Melles-Griot).  dc,  ac  and  phase  shift  signals  are  obtained 
using  the  standard  heterodyne  technique  controlled  by  Fast  Fourier  Transform  (FFT)  Labview 
routines.  Fluorescence  signals  are  obtained  through  an  700  ±  10  nm  interference  filter  placed 
in  front  of  the  detection  PMT.  For  each  experimental  configuration,  a  total  of  3x256 
measurements  (256  dc,  ac  and  phases  each)  can  be  made.  This  requires  about  8  minutes  at  the 
present  time. 

The  oxygen-sensitive  dye  used  (SCCN)  was  obtained  from  Frontier  Scientific  (Logan,  UT) 
and  used  without  further  purification.  The  peak  absorption  and  emission  for  this  dye  are  near 
635  and  694  nm,  respectively.  We  used  a  30  mm  diameter  cylindrical  solid  phantom  (1% 
Intralipid  +  10  mg/liter  SCCN  +  2%  Agar)  as  the  background  medium  which  had 

(Xa  =  0.005/mm  and  |la  =1.0/ mm .  A  glass  tube  (10  mm  outer  diameter,  9.5  mm  inner 

diameter)  was  inserted  as  a  target  at  different  positions,  as  shown  in  Figs,  la  and  lb.  The 
target  tube  contained  a  solution  with  the  same  amount  of  mixture  of  Intralipid  and  SCCN  as 
the  background.  Thus  both  the  target  and  background  had  the  same  scattering/absorption 
coefficients  and  fluorophore  concentration.  The  lifetime  contrast  was  introduced  by  varying 
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the  amount  of  oxygen  content  in  the;  target  and  background.  The  gas  delivery  system,  shown 
in  Fig.  lc,  was  a  simple  home-made  system  designed  to  obtain  only  three  differing  oxygen 
concentrations  where  the  oxygen  partial  pressure  (p02)  above  the  target  solution  was  0,  0.21, 
and  1.0  atm.  Since  we  found  that  lowering  the  oxygen  content  of  the  target  (0  atm)  gave  the 
greatest  lifetime  contrast  with  our  system,  we  used  0  atm  p02  above  the  target  and  1  atm  p02 
above  the  background  for  all  our  experiments.  In  these  experiments,  nitrogen  was  bubbled 
into  the  target  solution  for  a  half  hour  while  the  background  was  simply  open  to  air.  Finally, 
we  conducted  an  experiment  with  homogeneous  phantom  where  the  target  and  the 
background  had  the  same  oxygen  concentration  (1  atm).  This  experiment  was  aimed  to  see  if 
the  0.5  mm  thickness  of  the  target  tube  and  the  absence  of  agar  in  the  target  have  significant 
effect  on  our  results. 


Fig.  2.  (a)  Reconstructed  lifetime  image  with  lowered  oxygen  content  in  the  target  (3  o’clock). 

The  axes  (left  and  bottom)  illustrate  the  spatial  scale,  in  millimeters,  whereas  the  color  scale 
(right)  records  the  fluorescence  lifetime,  in  nanoseconds,  (b)  Lifetime  profile  along  a  horizontal 
cut  line  through  the  centers  of  both  the  target  and  background. 

3.  Results  and  Discussion 

The  finite  element  mesh  used  for  obtaining  the  results  shown  in  this  section  consisted  of  249 
nodes  and  448  triangle  elements.  The  final  images  reported  are  the  result  of  20  iterations. 
Figure  2a  presents  the  reconstructed  lifetime  image  when  the  target  was  positioned  at  3 
o’clock,  while  Fig.  2b  gives  the  lifetime  profile  along  a  cut-line  through  the  centers  of  both 
the  target  and  background  as  shown  in  Fig.  la.  As  can  be  seen,  the  target  is  clearly  detected. 
We  note  that  the  peak  lifetime  value  of  the  target  and  the  value  of  the  background  lifetime  are 
in  excellent  agreement  with  the  expected  values  of  1.14  and  0.80  ns  for  the  target  and 
background,  respectively  (the  expected  lifetime  values  of  the  dye  solution  at  0  and  1  atm  p02 
were  measured  using  a  standard  lifetime  spectrometer).  While  die  target  is  accurately  located, 
the  target  size  is  underestimated  in  this  case.  The  lifetime  image  and  its  profile  for  the  target  at 
a  deeper  location  (12  o’clock)  are  shown  in  Figs.  3a  and  3b.  While  we  note  that  the  target  is 
correctly  located  and  that  the  lifetime  values  of  both  the  target  and  background  are 
comparable  to  the  expected  values,  we  see  that  the  overall  quality  of  reconstruction  in  this 
case  is  worse  than  the  image  shown  in  Figs.  2a  and  2b.  We  believe  the  degradation  of  image 
quality  in  this  case  is  primarily  attributed  to  the  lower  signal-to-noise  ratio  (SNR)  available 
when  the  target  is  embedded  at  a  larger  depth  from  the  boundary.  This  low  SNR  case  also 
gives  more  artifacts  around  the  boundary  as  we  see  in  Fig.  3a.  The  recovered  lifetime  image 
for  the  homogeneous  case  is  depicted  in  Fig.  4  where  we  can  see  that  the  target  tube  and  agar 
concentration  used  do  not  affect  our  reconstructions  significantly. 

In  summary,  we  have  presented  reconstruction  of  lifetime  maps  based  an  oxygen-sensitive 
dye  in  turbid  media  using  ac  excitation  and  fluorescent  data.  We  have  demonstrated 
quantitative  imaging  ability  for  recovering  fluorescence  lifetime  in  centimeter  tissue-like 
phantoms.  Thus  oxygen  concentration  maps  can  be  obtained  for  biomedical  diagnosis.  It 
should  be  pointed  out  that  the  reconstruction  algorithm  presented  here  is  applicable  to  non- 
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singlet  delay  of  fluorescent  dyes  for  detection  of  low  oxygen  concentration  contrast  [4], 
Unlike  the  fluorescence  lifetime  imaging  (FLIM)  reported  thus  far  [19],  the  lifetime 
tomography  described  here  is  built  upon  the  model-based  reconstruction  method,  thus  suitable 
for  imaging  large  tissues  such  as  breast.  This  study  represents  our  initial  step  toward  in  vivo 
fluorescence  lifetime  tomography.  While  the  results  presented  are  promising,  we  realize  that 
improving  SNR  for  FLT  remains  as  a  challenge. 


Fig.  3.  (a)  Reconstructed  lifetime  image  with  lowered  oxygen  content  in  the  target  (12  o’clock). 
The  axes  (left  and  bottom)  illustrate  the  spatial  scale,  in  millimeters,  whereas  the  color  scale 
(right)  records  the  fluorescence  lifetime,  in  nanoseconds,  (b)  Lifetime  profile  along  a  horizontal 
cut  line  through  the  center  of  the  target. 


Fig.  4.  Reconstructed  lifetime  image  with  homogeneous  oxygen  content  in  the  target  and 
background.  The  axes  (left  and  bottom)  illustrate  the  spatial  scale,  in  millimeters,  whereas  the 
color  scale  (right)  records  the  fluorescence  lifetime,  in  nanoseconds. 
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